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Abstract 

Building on recent research into five-dimensional holographic models of QCD, we extend 
this work by including the strange quark with an SU{3)l x 5C/(3)i?, gauge symmetry in 
the five-dimensional theory. In addition we deform the naive AdS metric with a single 
parameter, thereby breaking the conformal symmetry at low energies. The vector and 
axial vector sectors are studied in detail and both the masses and decay constants are 
calculated with the additional parameters. It is shown that with a single extra degree 
of freedom, exceptional agreement with experimental results can be obtained in the light 
quark sector while the kaon sector of the vector and axial vector octet is found to give 
around 10% agreement with lattice results. We propose some simple extensions to this 
work to be taken up in future research. 
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1 Introduction 



Since 1973 QCD [1] has been the estabhshed theory of the strong interactions. Though un- 
derstood in the perturbative regime, there remains no analytic solution for the phenomenon of 
confinement for a general Yang-Mills theory. Effective chiral perturbation theory is constructed 
in order to describe the low energy interactions of matter fields. A systematic method to han- 
dle the non-perturbative regime of the strong interaction from first principles has been sought 
though without success so far. Inspired by the gauge/gravity duality [2, 3, 4] and the success 
of flavoured models from ten-dimensional supergravity [5, 6, 7, 8, 9], several holographic mod- 
els [10, 11, 12, 13, 14, 15, 16, 17, 18, 19] have been proposed recently with surprising success. 
The key point is that although QCD is not a conformal field theory and its string dual is still 
unknown, it approaches the conformal limit in the ultraviolet regime. Low energy properties 
such as confinement and chiral symmetry breaking are introduced in these models in a simple 
manner. The models discussed so far have all included quarks in SU (2) or SU (3) flavour groups, 
each with an identical mass. 

Since the mass differences between the mesons in the same octet are small compared with the 
scale Aqcd for the lowest vector and axial- vector mesons, it is interesting to study whether the 
same successful results of the isovector sector can be achieved after including the third flavour. 
We will show that indeed this is the case. In this letter we will extend the work of [10] to three 
flavours without assuming the flavour symmetry and calculate the mass spectrum and decay 
constants of l^Si vector mesons and l^Pi axial-vector mesons. We will not discuss the scalar 
sector and the U{1)a problem but will leave these questions for future research. 

The metric in the five-dimensional bulk with four-dimensional Poincare symmetry and a 
compact extra dimension is 

ds^ = 9{z)e{zjR - z)a^{z){dx''dx^ - dz^), (1) 

where 

1 

a^{z) = — (1 + dnZ^^) and 9{z) is the step function. (2) 

^ j=i 

The b*^ dimension corresponds to the energy scale in the four-dimensional theory. The metric 
has a double pole on the UV boundary {z = 0) and is asymptotically AdS^ in the z ^ limit. 
The IR boundary Zjr ~ 0(Ag|^^) is put in by hand in this model in order to introduce the 
mass gap and can be interpreted as the position of an IR D-brane. In this work we will consider 
both the simplest case, that is, the AdS bulk such that a{z) = ^ and the case where only ai is 
non-zero in the warp factor a{z). Since the AdS background is dual to A/" = 4 SYM, considering 
the deformed metric is natural for a model builder as QCD is not a conformal field theory. In 
the bulk the gauge group is chosen to be SU{3)l ® SU{3)ji, which is dual to the global chiral 
group on the UV boundary where QCD behaves approximately like a CFT. In addition to the 
gauge fields Lm{z,x) and Rm{z,x), a scalar field X transforming as {3l,3r) is also introduced 
to induce the chiral symmetry breaking. The non-renormalizable five-dimensional action is 
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where 



DmX — OmX — iLmX + iXRM, 

T T a j.a 

J^M — -^Af'^ , 

Lmn = OmLn - OnLm - i[LM,LN], similarly for i?, (4) 

where Tr[t",t''] = For our interest in this letter, we keep only terms up to quadratic 

order in fields. The effects of including higher order non-renormalizable terms is equivalent to 
deforming the metric from the AdS metric, which will also be considered in this work. However, 
the an terms are blind to flavours. The only source which violates SU{3) flavour symmetry comes 
from X. We deflne the vector and axial- vector gauge bosons to be Vm = ^{Lm + Rm) and Am = 
^{Lm — Rm) respectively. Following [10], we will choose the axial gauge = A^ = 0. According 
to the holographic recipe [4, 20], the non-normalizable mode of is dual to the conserved vector 
current in QCD while the normalizable modes give the flve- dimensional extension of the vector 
meson wavef unctions, similarly for A^ flelds. Also, the mass of the scalar fleld X is determined 
by the AdS/CFT correspondence. First, we consider the pure AdS case. After solving the 
equation of motion for the expectation value of X, the solution is 

{Xiz)) = ^iMz + J:z'), (5) 



where 



M = ,ta2<^ and j = 

.-0 z zj^ 



M is the quark mass matrix and S = {qq), which can be easily shown by matching the four- 
dimensional effective Lagrangian to the chiral Lagrangian [11]. S is responsible for the x^B 
in the chiral limit and we expect S oc 1 in this limit. That is, S ~ al + 0{M). Instead of 
considering possible boundary terms which determine the speciflc forms of M and S, we will 
parametrize the effects of these terms and choose M = diag{m^ m, rus) in the isospin limit for 
phenomenological interest and S = diag{c, c, Cg). With this choice, SU{3)l ® SU{3)ji is broken 
to SU{2) ® U{1). 

Substituting X = l^X)e'^'^^°"^°'^^'^^ back into the action, one can show that the flve-dimensional 
flelds and have block- diagonal z-dependent mass matrices My and Ma respectively where 

/ 03x3 

M^= \{m-ms + {c-Cs)z^f z^U^^ |, (6) 
\ 

and 

/ {m + cz'^fz^U.^^ \ 



Ml 



\{m + rUg + {c + Cs)z'^f z'^14^^4^ 
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yO \(2{ms + Csz'^f + {m + cz^fy z"^ J 

(7) 

V^''^'^, V^'^'^''^, and correspond to isovector, isodoublet, and isosinglet vector mesons in the 
octet of the quark model, similarly for the axial-mesons. Note that the mass of isovector and 
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isosinglet vector mesons, which are related to My, are the same and independent of {X). In 
the large Nc limit [21], where Zweig's rule can be shown to be exact, this result suggests the 
isosinglet in the octet to be the uj meson in the "ideal mixing" limit. The origin of the mass for 

is not from the Higgs mechanism but only from the explicitly flavour symmetry breaking 
term. In the '^'[/(S)^^^^^^^ limit, fields are massless. In this limit, all the vector mesons in 
the same octet have the same mass. 

The coupling constant is determined by matching the high energy expansion for the two- 
point function, calculated by the holographic recipe, to the same function calculated from the 
OPE in large Nc limit [22]. The result is [10] g'^ = ^j^- This result relates the coupling constant 
Nc in large Nc QCD to the one in its dual model. 

The vector meson, V, wavefunction, ^v, and the axial meson. A, wavefunction, which 
are the normalizable modes of four- dimensional Fourier transformed transverse field V^j_ and 
A'J^^ respectively, satisfy the following linearized equations of motion (no summation) 

[dl + (Ina(z)) + {q' - {glaizfM^U)] z) = 0, (8) 

and 

[d! + a, (Ina(z)) + {q' - {glaizfMlU)] z) = 0, (9) 

with boundary conditions dz^v{q, zm) = and $y(g, e) = 0, similarly for In the large 
Nc limit, the hadron masses can be obtained from the poles of the two-point function. This is 
equivalent to solving Eq.(8) and Eq.(9) with eigenvalues q"^ = my and q^ = m\ respectively. 
The decay constants F^a, which are defined by (0 | Jy A^i I ^) = ^"''^FyAe^ with the 
polarization vector, can be shown to be [10] 



where 



A^= / dza{z)\^v,Ai^)\ 
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The masses of pseudoscalar mesons P can be obtained by solving the following equations 
mixed by the normalizable modes of the fields n"-{x,z) and the longitudinal components of 
^4°, defined to be = ^^^"(Xjz), at q^ = mp with boundary conditions (p'^'^zm) = 

'^"''{zir) = (^°(0) = 7r°(0) = (no summation) 

{dl + dz (lna(^)) a,) r + gWiz) (Ml)^^ (tt^ - r) = 0, 

1 A - .ox 1 



{a'{z) [Ml + Ml)^^ d.n^) + a\z)q' (^[M^ + Ml) [n'^ - + [M^ - M]) 



aa 
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The decay constants of massless pseudoscalar mesons P° (which is a good approximation 
for pions) are given by [10] 

_ 1 a.A"(o, z) , 

fp^ - -J, U=., (12) 

with A^iQ, z) the solution of Eq.(9) satisfying ^'^'(O, zir) = and y4"(0, e) = 1. 



The results are shown in Table 2 and Table 3. It is well-known that without flavour symmetry, 
the l^Pi isodoublet state mixes with l^Pi isodoublet state. The mixing angle 6 is defined by 
parametrizing the physical state /ri(1273) and Jri(1402) 



\Ki{U02)) = \KiA)cose- \KiB)sm9, 
|Js:i(1273)) = \Kia) sine - \KiB)cose, 



(13) 



such that 

sin^^ 



^^-'i^l(1402) ^'^Xi(1273) 



In the numerical analysis we fit the value of the quark mass and condensate to the case of ideal 
mixing {6 = 45°). 

We now turn to the case where the warp factor is deformed from the AdS case. We consider 
the case that only ai is non-zero. In this case the solution {X{z)) is of the form {X{z)) = 
diag{Xu{z) , Xu{z) , Xs{z)) where Xu{z) and Xs{z) approach their AdS limit respectively when 
2; — *■ 0. All the qualitative discussions above still hold and our numerical results are shown 
in Table 3 and Table 4. ai is expected to be of order 0{g1m?^ from matching the two-point 
function to the same function in the large A^^^ QCD. Our numerical result for ai agrees with 
this. 



2 Numerical Analysis 

Precisely as in [10] the value of the IR cutoff is calculated using the equation of motion of the p 
meson, which is independent of the quark mass. The value of this cutoff is however a function 
of the deformation parameter, «i, which is tuned such that the p decay constant is closer to the 
experimental value. Note that by tuning ai, it is possible to obtain exactly the correct value for 
this decay constant. In this work a semi-global fit has been used so that the value of the decay 
constant is good but not exact, such that the predictions of other observables are closer to the 
experimental values. This method is similar to the global fit used in [10] but here, because we 
have three more free parameters, a full global fit is unmanageable. This presumably means that 
a full survey of the parameter space would give even better results than those presented here if 
such a fit were computable. 

Having calculated the IR cutoff, the excited states of the p meson are found by studying the 
higher modes of the appropriate five- dimensional field. 

In [10], a perturbative expansion is used to prove the Gell-Mann-Oakes-Renner relation [23]. 
We show here that an alternative numerical approach can be used to calculate the values of 
the quark mass m and condensate c. The value of the pion decay constant is a function of 
the quark mass and condensate. Similarly the boundary behaviour of the five- dimensional field 
corresponding to the pion is a function of these parameters (with set to in the pseudoscalar 
equation of motion). We can plot both of these quantities as a function of the quark mass and 
condensate. Each plot should have a line in (m, c) space corresponding to the physical values 
(the correct value of and the correct boundary behaviour respectively). The point at which 
these lines intersect therefore gives the correct value of m and c to provide the physical pion 
mass and pion decay constant. These "physical contours" are shown in figure 1 to illustrate 
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that this method produces almost identical values to the perturbative approach [10]. In both 
the undeformed and deformed AdS cases it can be shown that the GMOR relation holds to 
order as expected. 




Figure 1: Intersection of the physical contours for 771,7^ and /^r in the (m,c3) plane. 

The mass of the ai is calculated by solving the axial vector equation of motion with the 
appropriate boundary conditions. It is found that in the case of deformed AdS, the mass 
and decay constant values for this axial vector particle are in remarkably close agreement with 
experiment. 

For the kaon sector, we perform a global fit in the space of and Cg by studying the mass 
of the vector and axial vector Ji'-mesons corresponding to the i^*(892) and Ki^ mesons. It 
is possible to get within a few percent of the experimental value (in the case of ideal mixing) 
though it may be possible by tuning the deformation parameter and the values of m and c 
correctly to get these masses exact. 

This concludes the calculation of the free parameters of the theory. Using these values, the 
decay constants for the isodoublet and isosinglet sectors are calculated using Eq.(lO). 

We have shown by numerical analysis that even without a full global fit, the vector sector 
is in close agreement with the experimental results. Compared with the results of the QCD 
sum rules, this holographic action appears to model the vector and axial vector spectrum of 
the strong force to high precision. A global fit would almost certainly reduce the average error 
to the order of a percent or two, however as explained above such a calculation would require 
considerable computing time. 

We provide in Tables 1, 2 and 3 the values calculated in our analysis. 
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Parameter 


Value (MeV) with ai = 225^ 


Value (MeV) with ai = 


1 -ft 


335.6 
2.012 


322.6 
2.29 


1 

rris 


346.25 
90 


327.25 
100 


1 

„3 


457.53 


232.5 



Table 1: Values chosen for the free parameters for both pure AdS^ and the deformed geometry. 
Note that in the case of ai 7^ 0, we have performed a semi-global fit over pairs of variables. 



Observables 


Value (MeV) with ai = 225^ (% error) 


Value (MeV) with ai = (% error) 




139.6* 


139.6* 


fn 


92.4* 


92.4* 




1238 (0.24) 


1363 (10.4) 




434.92 (0.44) 


486 (12.2) 




1339* 


1339* 




436 (v/i"/<i(i400) ~ 454t) 


450 (v/i^i^i(i4oo) ~ 454t) 




1389 


1339 


\/Fa3 


430 


443 



Table 2: Axial sector results for both pure AdS^ and the deformed geometry using the param- 
eters in table 1. Note that the ai = results have been calculated using the same numerical 
techniques as in [10]. The case of ideal mixing {6 = 45°) has been chosen to fix the strange 
quark mass and condensate values. Experimental values are chosen as the midpoint of those in 
[24]. The decay constant of the ai is compared with the lattice result [25]. * indicates that this 
value is used to fix the free parameters, all other values are predictions, f results are taken from 
[26]. The axial vector meson A^ corresponds to the isosinglet meson in the octet. 



Observables 


Value (MeV) with ai = 225^ (% error) 


Value (MeV) with ai = (% error) 




775.8* 


775.8* 




335.5 (2.8) 


329 (4.5) 




1830 


1781 




635.38 


616.4 


niK* 


793 (11)* 


799 (10.4)* 




330.8 (11^) 


329 (lit) 




mp 


nip 









Table 3: Vector sector results for both pure AdS^ and the deformed geometry. The p' is the first 
excited state of the p meson. Experimental values are chosen as the midpoint of those in [24]. 
* indicates that this value is used to fix the free parameters, all other values are predictions, f 
results are taken from lattice predictions [27] though these have large uncertainties. The vector 
meson V3 corresponds to the isosinglet meson in the octet. 
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3 Discussion and Conclusion 



We have extended the work of Erhch et al [10] to include an SU{3) flavour group with a distinct 
strange quark mass. We also study the effects of a deformation of the metric by the simplest 
possible factor. This deformation is sufficient to provide a non-conformal region in the effective 
four-dimensional theory. We show that the effect of adding this extra parameter is enough to 
give very accurate (e^ms < 2%) masses and decay constants for the light quark sector. 

We extend the formalism with the addition of a strange quark and calculate the vector and 
axial vector masses and decay constants in the s quark sector. The values of the strange quark 
mass and condensate are tuned to give the Ki^ and K* masses to within a few percent of the 
experimental values (though we assume ideal mixing). This fit is a function of the deformation 
parameter, cti, so that with a full global fit it may be possible to lower the average error, which 
in this analysis is ~ 7%. With the mass and condensate values fixed, we calculate the decay 
constants which are compared with lattice data providing a good agreement. 

Though in the sectors described in this paper an excellent agreement with experiment and 
lattice is obtained, it is found that the strange pseudoscalar (the kaon) gives a surprisingly poor 
fit for the parameter values chosen here. It is also possible to tune m<j and Cg such that the kaon 
mass is correct giving a good fit for the K* but a poor result for the Ki. 

Note that the values of m„ and c have been fixed without reference to the kaon sector. It 
may well be possible to tune these values such that the kaon sector provides a better fit to the 
data. In this letter we merely wish to illustrate that these extensions to the work of Erlich et al 
[10] can provide reasonable fits in the strange quark sector. We leave a full global fit for further 
calculation. 

There are several possible extensions to this work. It would be interesting to study the 
mass splitting in the isovector spectrum by the addition of different mass u and d quarks and a 
four-dimensional U{1) gauge interaction. It would also be interesting to study the addition of 
an operator corresponding to baryons in these models. Several other interesting directions are 
noted in [10] which should guide us even closer to a true holographic dual of QCD. 
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